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Abstract. We present the latest results on the search for bottom baryon states St using 
■ ~ lfb -1 of CDF data. The study is performed with the world's largest sample of fully 

C ' reconstructed A® decays collected by CDF II detector at y/s = 1.96 TeV in the hadronic trigger 

path. We observe 4 new states consistent with X'£*- )± bottom baryons. 

X 

. 1. Introduction 

High energy particle colliders provide a wealth of experimental data on bottom mesons. However, 
only one bottom baryon, the A®, has been directly observed [U [31 0] . 

Heavy baryons containing one heavy quark and a light diquark became a nice 3-body 
laboratory to test QCD models. In the limit of heavy quark mass txiq — * oo, heavy baryons' 
properties are governed by the dynamics of the light diquark in a gluon field created by the 
heavy quark acting as a static source. The heavy baryon like A® can be considered as a "helium 
atom" of QCD. In this Heavy Quark Symmetry (HQS) approach [5] at the heavy quark limit a 
heavy quark spin does not interact with the gluon field, the spin decouples from the degrees of 
freedom of the light quark and the quantum numbers of the heavy quark and the light diquark 
are separately conserved by the strong interaction. Consequently the light diquark momentum 
jqq = s qq + L qq , the heavy quark spin sq and total momentum JQ qq = sq + j qq are considered 
as good quantum numbers. Based on the HQS principles an effective field theory was constructed 
where corrections can be systematically included in the perturbative expansions. The theory 



X 



was named as Heavy Quark Effective Theory (HQET) (see [6] and references therein). 

For the bottom Q = b baryons with a single heavy quark and two light ones (see Table [1]) 
the bottom quark spin, sq = is combined with the light diquark momentum j qq comprised 
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of spin Sqq — 0+ © 1+ and its angular momentum Lqq- The baryons with a diquark having 
s qq = + and isospin 1 = are called A- type, while the states with s qq = 1 + and isospin 1=1 
are called E- type. The ground state A® baryon has I = 0, J p = \ .A doublet of ground E- 

like bottom baryons comprises E b with I = 1, J p = | + and E b with I = 1, J p = | + . 

The combination of an orbital momentum L qq = 1 of a diquark with its spin of s qq = 0+ 
adds to the spectroscopy a number of excited P- wave bottom A -states. The lowest lying orbital 
excitations are Af with J p = \~ ', |" [3E1E]. 

Table 1. Bottom baryon A- and E- states and their quantum numbers. The \q\q2\ denotes a 
pair antisymmetric in flavor and spin. The {(71(72} denotes a pair symmetric in flavor and spin. 
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Theoretical expectations for ground bottom baryon states are summarized in Table The 
calculations have been done with non-relativistic and relativistic potential quark models, 1/N C 
expansion, quark models in the HQET approximation, sum rules, and finally with lattice QCD 
models [TP] . 

In a physics reality hiq is finite and a degeneration of a {E b , E^} doublet is resolved by 
a hyperfine mass splitting between its states. There is also an isospin mass splitting (see 
Table [2]) between members of E^ and E£ isotriplets [HI [121 El El- As it was pointed 
out |14j the value of the isospin splitting within E^ triplet does differ from E^ triplet, namely 
- El~) - (E£ - J7 6 ~) = 0.40 ± 0.07MeV/c 2 . This number contributes to systematic 
uncertainty of our experimental results (see Section [5]). 



Table 2. Mass and width predictions for E^ 
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see below 





According to HQS the physics of pion transitions between heavy baryons is governed by the 
light diquark. The one or two pions are emitted from the light diquark while the heavy quark 



propagates unaffected by the pion emission process. Various pion transitions of bottom baryons 
into the lower ground states are summarized at the Figure [H The mass predictions for the S- 
wave (i.e. ground state) E- like baryons show that there is enough phase space for the both Eb 
and E£ to decay into A® via single-pion emission. The two excited (P-wave) states might 
decay into A® via two-pion transitions provided a sufficient phase space. 
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Figure 1. Transitions of ground 
states Ef,, E£ into a low lying 
ground state A® via a single-pion 
emission in a P-wave as the diquark 
with j qq = 1 + is converted to one 
with j qq = 0+. If Af states 
have masses sufficiently higher than 
the decay thresholds, the two-pion 
transitions into A® or a single-pion 
decays into Ef, or El states are 
possible (see also a discussion [7, 
El |9j). Al° single-pion transitions 
to A® are forbidden due to isospin 
conservation and also by parity (for 
the higher J p = | state). 



It is important for our experimental expectations to understand the natural width of E^ 

baryons. As we expect that Zv masses lie within (180 - 210) MeV/c 2 above A° b and well above 
a threshold for a single-pion mode, we would expect that the single-pion P-wave transition will 
dominate the total width [I5j [16] . The authors [16] find 

1 M Aq 2-3 

r^..^ = — ^— |/ p | \vA , 

where p n is the three- momentum of soft K£q, f P = 9aI ' fin In = 92MeV and qa is the axial 
vector coupling of the constituent quark for the nucleon. A fit of this formula to the known 
PDG width measurements [T7] of charm states E c and E* yields qa = 0.75 ± 0.05 which is in 
excellent agreement with qa = 0.75 numerical theoretical value for the nucleon. Using the fit 
results we have estimated r(Z^), see the Table [2j The error of the fitted qa contributes as a 
systematic uncertainty to our experimental measurements (see Section [5]). 



2. Principle of the Analysis 

The topology of the event with E b state produced in Tevatron collisions is demonstrated in 
Figure EJ The E b candidates are searched in the decay chair0: 

• Strong decay S^*^ — > A^tt^ with both A® — > A+ir~ and its daughter A+ — > pivT~7r + in 
weak decay modes. 

To remove a contribution due to a mass resolution of each A® candidate and to avoid absolute 

mass scale systematic uncertainties, the U b candidates are reconstructed in the mass difference 
Q-value spectra defined as 

Q = M(A° b Tr± b ) - M(A° b ) - Mpdg^). 

The narrow signatures are searched for in the Q-spectrum constructed separately for every charge 
state of JC^ ^ candidates. The subsample of S b contains A®ir~ and A b ir + combinations from 
the decays of the particles S b *^ and the antiparticles , respectively. The subsample of 

U b *^ + contains ^5vr + and A®ir~ combinations from the decays of the particles ^* and the 
antiparticles U b , respectively. 




Figure 2. Sketch of the event topology 
of a U), produced in the CDF detector [18] . 
The bottom baryon U b is produced at the 
collision origin, i.e. primary vertex of the 
event. The Ub decays strongly into A® with 
soft 7r^ emitted. Due to a fast nature of the 

^b 

strong decay the soft pion originates from the 
primary vertex. The A® and its daughter A£ 
decay weakly with a secondary decay vertex 
measured in SVX II [19]. The tracks with 
a common secondary vertex are displaced 
relative to a primary vertex and the proton 
from A~£ decay and ir~ from A® decay most 
likely contribute to the CDF hadron Two 
Track Trigger (see Section [3] and |20j). 



The S b *^ signal region at Q- value spectrum is defined as 30MeV/c 2 < Q < 100MeV/c 2 , 

based on the theoretical expectation (see Table [2]). We pursued a blind analysis and developed 

(*)± 

the E b selection criteria using only the pure background sample in the upper and lower 
sideband regions of OMeV/c 2 < Q < 30MeV/c 2 and 100MeV/c 2 < Q < 500MeV/c 2 . The signal 
was modeled by a PYTHIA [H] Monte Carlo. 

3. Triggers and Datasets 

Our results are based on data collected with the CDF II detector [18] and corresponding to an 
integrated luminosity of ~ l.lfb~ . As pp collisions at 1.96 TeV have an enormous inelastic 
total cross-section of ~ 60 mb, while b- hadron events comprise only ~ 20 /zb (\rj\ < 1.0), 
triggers selecting b- hadron events are of vital importance. Our analysis is based on a data 

2 Unless otherwise stated all references to the specific charge combination imply the charge conjugate combination 
as well. 



sample collected by a three- level Two displaced Track Trigger. It reconstructs a pair of high 
Pt > 2.0GeV/c tracks at Level 1 with the CDF central tracker and enables secondary vertex 
selection at Level 2. This requires each of the tracks to have an impact parameter measured 
by the CDF silicon detector SVX II [19} [20] to be larger than 120 //m. The excellent impact 
parameter resolution of SVX II makes this challenging task possible. The trigger proceeds with 
a full event reconstruction at Level 3. The Two displaced Track Trigger is efficient for heavy 
quark hadron decay modes (see Figure E] and its caption) . 

4. Event Selection 

The candidates of the basic state in our analysis, A®, have been reconstructed in the mode 
A b — > A^n~ with A£ — ► pK~ir + . The Two displaced Track Trigger requirements (see 
Section [3j) are confirmed offline for each A® candidate. The Charm Af and bottom A® 
candidates have both been subjected to 3-dimensional vertex fits. The collection of the fitted A~£ 
candidates has been confined to a mass range of m p ^ G ± 16 MeV/c 2 [17J . To suppress a prompt 

background we apply a cut on the proper decay time cr(A b ) > 250 /xm with its significance 
cr(A b ) I 'a CT > 10. The proper decay time of A£ with respect to the A® vertex is required to be 
-70 urn < ct(A+ «- A° b ) < 200 fim. We define the topological quantities as cr = L xy ■ tuaqIvt 
and L xy = D xy ■ pr/pr where D xy , p~r are the vectors of corresponding distance or momentum 
in a transverse plane. To reduce combinatorial background and contribution from partially 
reconstructed modes the impact parameter do(A b ) is also restricted to be below of 80 /urn, where 
do = \D xy x p~t/pt\- The kinematic cuts for A® and A+ candidates, pr(A®) > 6.0GeV/c and 
Pt(A+) > 4.5GeV/c are applied as well. 

The powerful A® — > A^ir~ signal is shown in Figure [3] with a binned maximum likelihood fit 
superimposed. The background from physical states contributing to the left side of the signal is 
analyzed with Monte Carlo simulations. The fit to the invariant A+tt~ mass distribution yields 
3125 ± 62(stat) of A° b -> yl+vr" candidates. We posses the world's largest A® sample. 
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Figure 3. Fit to the invariant 
mass of A® — > A+ir~ candidates. 
Fully reconstructed A® decays such 
as A° b -> A+7T- and A° b -> A+K~ 
are not indicated on the figure. The 
A b signal region, 5.565 GeV/c 2 < 
m(A+iT-) < 5.670 GeV/c 2 , consists 
primarily of A b baryons, with some 
contamination from B mesons and 
combinatorial events. The left 
side band is enriched by partially 
reconstructed A b decays like A b — ► 
yl* + 7r~ and by fully reconstructed 
4-prong B- meson decays like B° — > 
D+ir-, D+ K-ir + ir + . The 

right side band consists from a 
combinatorial background. 



5. Sb Candidates and Signals 

Following a method outlined in a Section [2] the A b candidates (see a Section H]) from a signal 
region of [5.565, 5.670] GeV/ c 2 have been coupled with the pion tracks 7r| to create 



candidates and the pairs A^n^ have been subjected to a common vertex fit. 

(*) 

The analysis cuts are optimized according to a blind analysis technique using the upper 
and lower sideband regions (see a Section d]) of experimental Q- value spectrum while the signal 
is modeled by a PYTHIA |21j . The following kinematic and topological variables are used : 
PT(^b)y the soft pion track impact parameter significance | <^o / I (^a)) an d the polar angle of 
the soft pion in a U^- rest frame, cos 9* (iV£ b ) = ps b ■ Pt/(\Ps b \ ' KD- A figure of merit is defined 
as €(Smc)/VB, where e(S , Mc) is the signal efficiency measured in the Monte Carlo sample and 
B is the background in the signal region estimated from the upper and lower sidebands. The 
maximum of the figure of merit is reached for pxi^b) > 9.5GeV/c, \do/ad \ {^s h ) < 3.0, and 
cos 9* > -0.35. 




Q = m(A?7r) - m(A?) - n\ (GeV/c 2 ) Q = m(A°7i) - m(A°) - (GeV/c 2 ) 

a) b) 

Figure 4. a) The three background sources described in the text and their sum are shown 
superimposed on the Q distributions, with the signal region blinded. The top plot shows the 
A®7t£ distribution, which contains all . The bottom plot shows the A®ir~£ distribution, with 

all . The sideband regions are parameterized with a power law multiplied by an exponential. 
The percentage of each background component in the A® signal region is computed from the A® 
mass fit (see Figure E|, and is: N(A° b ) 90.1%, N(B) 6.3% and JV( combinatorial) 3.6%; 
b) The excess in the A b TT^ b (top plot) subsample is 148 candidates over 268 expected background 
candidates, whereas in the A®tt~£ (bottom plot) subsample the excess is 108 over 298 expected 
background candidates. 

The Q- value spectra with blinded signal region are shown in Figure 0^, with its detailed 
caption. In the search, the dominant background is from the combination of prompt A® 
baryons with extra tracks produced in the hadronization process. The remaining backgrounds 
are from the combination of hadronization tracks with B mesons reconstructed as A® baryons, 
and from combinatorial background events. 

Upon unblinding the Q signal region in both spectra we observe an excess of events over the 
background as shown in Figure Hb with the details explained in the caption. 



Next we perform a simultaneous unbinned maximum likelihood fit to the A^ir^ 



subsamples for a signal from each expected E^ state plus the background, referred to 
as the "four signal hypothesis." Each signal consists of a Breit-Wigner distribution convo- 
luted with two Gaussian distributions describing the detector resolution, with a dominant 
narrow core and a small broad component for the tails. The natural width of each Breit- 
Wigner distribution is computed from the central Q- value (see a Section [1] and |16j). The 
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fit shown in Figure 5 results in the yields N 
and 
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Table 3. Likelihood ratios calculated for the 
alternative signal hypothesis with respect to 

the one of four E^ states. The strength of 
the signal hypothesis is further given by the 
likelihood ratio, LR = L/L a \ t , where L is the 
likelihood of the four signal hypothesis and 
L a it is the likelihood of an alternate 
hypothesis. 
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Figure 5. Simultaneous fit to the A^tt^ 



The alternative hypothesis is estimated using 
all systematic variations of the background 
and signal functions. The variation with the 
largest value of L a \ t corresponding to the least 
favorable hypothesis is taken. 



(top) and A®tt~£ (bottom) spectra for E^*^ 
candidates shown on a range of 
Q S [0, 200] MeV/c 2 . 

Systematic uncertainties on the mass difference and yield measurements fall into three 
categories: mass scale, E^ background model, and E^ signal parameterization. The mass 
scale is determined from the difference in the mean of the narrow resonances m(D* + ) — m(D°), 
m(E®' ++ ) — m(A^), m(yl c (2625) + ) — m(A^) between data and PDGfTTj, The uncertainties on 
the background come from the assumption on the sample composition of the A® signal region, 
the normalization and functional form of the A® hadronization background. The systematic 

effects related to assumptions made on the E^ signal parameterization are: underestimation 
of the detector resolution in Monte Carlo, the accuracy of the natural width prediction 
from [16], and the fit constraint that {Q E *+ — Q E +) = (Q s *- — Q^-) [14]. All systematic 

b b b b 

uncertainties on the Q „(*)± mass difference measurements are small compared to the statistical 
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uncertainties. The final results [22] for the signal yields, including systematic errors, 
33l^ (statOig (syst.), N E - = 62+^ (stat.)^ (syst.), A^*+ - Q ^+ 17 '^^ 
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The final results [22] for the masses are Q, 



48.2±a;jj (stat.)lSi (syst.) MeV/c 2 , Q E - = 55.9lo[J (stat.)loi ( s y st M eV/c 2 , and A s * = 

b " 

21.5±f g (statOloj (syst.) MeV/c 2 . Using the CDF II measurement of m A o = 5619.7 ± 
1.2 (stat.) ± 1.2 (syst.) MeV/c 2 [23], the masses of the four states are [22]: 

mr+ = 5807.5±2:1 (stat.) ± 1.7 (syst.) MeV/c 2 , m r = 5815.2±J;g (stat.) ± 1.7 (syst.) MeV/c 2 , 
m s ,+ = 5829.01^ (stat.) ±1.7 (syst.) MeV/c 2 , m«- = 5836.7l 2 j (stat.) ±1.7 (syst.) MeV/c 2 , 

where the systematic uncertainties are now dominated by the total A® mass uncertainty. 

The significance of the signal is evaluated with two methods: using statistical Monte-Carlo 
pseudo-experiments and comparing the likelihoods of the default four signal hypothesis with 
pessimistic alternate ones. The randomly generated background samples are fit with the four 
signal hypothesis. The probability for background to produce the observed experimental number 
of signal events or more is found to be less than 8.5 x 10 -8 , corresponding to a significance of 
greater than 5.2 a. The results on study of likelihood ratios are summarized in a Table 3 and 
its detailed caption. 

6. Conclusions 

In summary, using a sample of ~ 3100 A\ — > A~£ir~ candidates reconstructed in l.lfb -1 of 
CDF II data, we search for resonant ^vr 1 * 1 states. We observe a significant signal of four states 

(*) 

whose masses and widths are consistent with those expected for the lowest-lying charged 
baryons: , , and . This result represents the first observation of the Zli baryons. 
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